Breathing Spots in a Reaction-Diffusion System 
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A quasi- 2- dimensional stationary spot in a disk-shaped chemical reactor is observed to bifurcate to 
an oscillating spot when a control parameter is increased beyond a critical value. Further increase 
of the control parameter leads to the collapse and disappearance of the spot. Analysis of a bistable 
activator-inhibitor model indicates that the observed behavior is a consequence of interaction of the 
front with the boundary near a parity breaking front bifurcation. 

PACS numbers: 47.54.+r, 82.20.Mj, 82.40. Ck 

Oscillations in spatially extended chemical systems are 
often the result of underlying oscillating dynamics of the 
local chemical kinetics jjj. In systems with nonuniform 
spatial structures, oscillations may also be driven by dif- 
fusion. Spiral waves and breathing spots in excitable 
and bistable media are examples of oscillatory behav- 
iors where the local kinetics without diffusion converge 
to stationary uniform states, while spatial structures un- 
dergo oscillations. Chemical spirals have been observed 
and studied for more than two decades ||] , but breathing 
spots have not been previously observed although they 
have been found in numerical and analytical studies of 
activator-inhibitor models ||||]. 

Figure 1 shows an example of the breathing spots ob- 
served in our study of a ferrocyanide-iodate-sulfite (FIS) 
reaction || in a quasi-2-dimensional reactor. The breath- 
ing motion arises as a control parameter is increased and 
an initially stable circular front (the spot boundary) be- 
comes unstable. Further increase in control parameter 
eventually leads to the front rebounding from the cell 
boundary and propagating inward until the spot collapses 
and disappears. The breathing motion is interpreted as 
transitions between left and right propagating fronts near 
a parity breaking front bifurcation |(J . The rebound phe- 
nomenon leading to spot collapse is attributed to cross- 
ing the front bifurcation as the control parameter is in- 
creased. We will first describe the experimental system 
and then present the observations and the interpretation 
of the results in terms of a model reaction-diffusion sys- 
tem. 

The chemical patterns form in a thin gel layer that 
allows reaction and diffusion processes but prevents con- 
vection. The apparatus is similar to that used by Lee et 
al. [Q. A polyacrylamide gel layer (25 mm diameter, 0.3 
mm thick) is in contact with a well-stirred reservoir (2.8 
ml volume) that is continuously fed with reagents of the 
FIS reaction. Reagents are fed first to a premixer (1.0 
ml volume) in two streams, one with H2SO4 and NaI03 
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FIG. 1. An oscillating circular spot at (a) minimal size 
and (b) maximal size, (c) Time evolution of the spot diame- 
ter. The input concentrations to the reservoir in contact with 
the gel reactor are [H 2 S0 4 ] = 3.35 mM, [IO~] = 75.0 ml, 
[SO^ 2 ] = 89 mM, and [K 4 Fe(CN) 6 • 3H 2 0] = 20 mM. The 
flow rate is 179 ml/h. 

and the other with Na 2 S0 3 and K 4 Fe(CN) 6 • 3H 2 0. The 
output of the premixer is fed to a stirred reservoir that 
is in contact with the gel layer. The reservoir diameter 
is 22 mm; thus the outer 1.5 mm width edge of the gel 
is not in contact with the reservoir. The entire system 
is immersed in a water bath maintained at T = 30° C. 
The side of the gel opposite to the chemical reservoir 
is a window through which the gel is illuminated with 
blue light (400-440 nm), and the patterns are viewed us- 
ing a video camera. The system is studied in a range 
in which the homogeneous FIS reaction has two stable 
states, one with low pH (about 4) and the other with 
high pH (about 7) (^J?]]. In the observed patterns the 
black and white regions correspond to the low and high 
pH states, respectively [Q. 

We will now describe experiments that yield oscillat- 
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ing spots. At low flow rates the gel reactor is uniformly 
black; at high flow rates, uniformly white Q. When the 
system is switched on at intermediate flow rates, a black 
spot emerges in the center of the reactor. The spot is ini- 
tially irregular but evolves to an circular spot centered in 
the reactor. Above a critical flow rate (about 150 ml/h 
for the reagent concentrations used; see Fig. 1), the circu- 
lar spot oscillates in size as it approaches its asymptotic 
stationary state. Beyond a higher critical value of the 
flow rate (158 ml/h), the circular spot becomes unstable 
and begins to oscillate in size, as Fig. 2 illustrates. The 
oscillations are nearly sinusoidal just beyond the onset 
of instability (Fig. 2(a)), while well beyond onset they 
become relaxational (Fig. 2(b)). Measurements of the 
amplitude A of the oscillation as a function of flow rate 
indicate that the transition is a Hopf bifurcation: A 2 in- 
creases linearly with distance above transition, as shown 
in Fig. 2(c). 
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FIG. 2. Time evolution of the cross section of a spot, (a) 
Sinusoidal oscillations just beyond the onset of instability 
(F — 160 ml/h). (b) Relaxation oscillations far beyond the 
instability onset (F = 200 ml/h). (c) The square of the spot 
oscillation amplitude as a function of flow rate. 



Beyond a yet larger flow rate (260 ml/h), a shrinking 
spot does not stop shrinking at a minimum size but in- 
stead continues to shrink until the spot disappears. As 
the spot collapses, a new black ring emerges near the 
outer edge of the reactor, creating two new fronts, as 
Fig. 3(d) illustrates. The inner front travels inward while 
the outer front initially travels towards the boundary but 
then rebounds and travels inward. The ring then col- 
lapses to a black spot that shrinks and disappears. The 
space-time diagram in Fig. 3(d) shows that this whole 
process is periodic. 
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FIG. 3. Periodic emergence of black rings near the bound- 
ary leading to spot collapse near the center of the reactor. 
(a)-(c) The dynamics of a single period, viewed at 80mm, 
85mm, and 93mm, respectively; the time origin is arbitrary, 
(d) The space-time evolution of a cross section of the reactor 
image. The flow rate is 280 ml/h. 

We interpret these observations as interactions of the 
chemical front with the reactor boundary in the vicinity 
of a parity breaking front bifurcation. The parity bro- 
ken front states correspond to a black state invading a 
white state (black-white front) and a white state invad- 
ing a black state (a white-black front). Although the 
two fronts connect the same states, they differ in their 
inner structures and consequently in their direction of 
propagation. To develop this interpretation we consider 
a model of a bistable reaction-diffusion system that ex- 
hibits pattern phenomenology similar to that observed in 
the FIS reaction ||, and is simpler for analysis than the 
G asp ar-S ho waiter model for this reaction Q. The model 
equations in one space dimension are 

u t = u - u 3 - v + u xx , (la) 
v t = e (u - aiv - a ) + Sv xx , (lb) 

where the subscripts x and t denote partial derivatives. 
For a fixed a± > 1 there is a parameter range where the 
system is bistable; it has two coexisting stable stationary 
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uniform states and (u-,V-). The parameter ao 

will be associated with the flow rate F, and the two states 
v+) and (u-,V-) with the white (high pH) and black 
(low pH) states, respectively. 

For ao = and fixed S the system (|l]) exhibits a pitch- 
fork front bifurcation (also known as the nonequilibrium 
Ising-Bloch (NIB) bifurcation) as the parameter e is de- 
creased past a critical value, e c . For e > e c there is a 
single, stationary front solution connecting (u+ : v+) at 
x — — oo to V-) at x = oo. At e = e c the stationary 
front solution becomes unstable and a pair of counter- 
propagating front solutions with velocities c oc ±^/e c — e 
appear. These are the parity broken front states corre- 
sponding to the black-white (c < 0) and the white-black 
(c > 0) fronts. When ao ^ the pitchfork bifurcation 
becomes imperfect, i.e., unfolds into a saddle node bifur- 
cation where, at e = e c (ao), a stable-unstable pair of front 
solutions appears in addition to the stable front solution 
that already exists. 

The front bifurcation can also be traversed by varying 
other parameters, in particular by increasing ao < 0. We 
will investigate the effect of a no-flux boundary on the 
dynamics of a front as ao is increased. Using a singular 
perturbation approach with e/S as a small parameter, we 
derive a relation between the front velocity, c, and the 
distance from the front to the boundary, d. The Hopf 
bifurcation observed in the experiment will be associated 
with the c — d relation becoming multivalued. 

In a moving coordinate system (1) becomes 



u zz + cu z + u — u - v = , 
5v zz + cv z + e (u — a\v — ao) = , 



(2a) 
(2b) 



where Xf(t) is the position of the narrow front structure, 
(t), and c = Xf is the front velocity. The 
boundary conditions are (u,v) — > (u-,v-) as z — > oo 
and (u z ,v z ) = (0,0) at z = —a 7 , where d is the distance 
from the front to the boundary. In obtaining ^) we as- 
sume that the front velocity Xf = c is small so that the 
explicit time dependence in the moving frame can be ne- 
glected. The front velocity can be controlled by varying 
ai; increasing a\ leads to lower velocities. 

We first solve equations ([|) in the front, or "inner", 
region. Letting /i = e/5— >0at finite rj = y/eS, we obtain 
the equation u zz + cu z + u — u 3 — Vf = 0, subject to 
the boundary conditions u — > u+(vf) as z — > — oo and 
u — > U-(vf) as 2 — > oo. Here, is the value of u at the 
front, and ^±(^/) are the largest and smallest roots of 
u — u 3 — Vf = 0. Solving the inner problem yields 

V2 

v f = — — c . 



(3) 



Another relation between c and Vf is obtained by solv- 
ing the equations in the regions to the left and to the 
right of the front, the "outer" regions. Rescaling the co- 
ordinate system according to £ = ^JJiz and letting \i — > 
gives Q 



g (v - v+) 

2/ 



^cc ^ — V C — r i (v — V-) = 



C<o, 
C>o, 

v/, v(oo) 



with the boundary conditions, v(0) 
and VQ(—yJJid) = 0. Here, a 2 = a\ + 1/2 and v± = 
(±1 — ao)/g 2 . Solving this boundary value problem and 
matching the derivatives of i? at £ = we find 
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(4) 



where a = \/ c 2 + 4eJg 2 . 

Equating (^) and (||) gives an implicit relation between 
c and d. In the limit d — > oo this relation reproduces 
the front bifurcation line j|, which we may write as 
&o = a>ob(e)- Fig. 4(a) shows a graph of the c — d re- 
lation far into the single front regime; Fig. 4(b), near 
the front bifurcation; and Fig. 4(c), beyond the bifurca- 
tion. The figures also show trajectories representing the 
front dynamics as obtained by direct numerical solution 
of equations (|l|). 

The monotonic velocity relation and the corresponding 
trajectory in Fig. 4(a) describe the approach of a black- 
white front from the far right (large x) to the boundary at 
x = Xf — d. The solution converges to a stationary front 
at some d = do- The negative slope at (c, d) = (0,do) 
implies stability of the stationary solution. We associate 
with this scenario the formation of a black spot observed 
at low flow rates in the experiment. 

As ao is increased, the slope of the c — d relation at 
(c, d) = (0, do) increases in absolute value and at a crit- 
ical point, aoc = — 1 + 2g 3 V2eJ/3, diverges to infinity. 
Beyond ao c the slope is positive and the c — d relation is 
multivalued in some range of distances from the bound- 
ary (Fig. 4(6)). The critical point ao c corresponds to the 
onset of oscillatory front motion, and explains the Hopf 
bifurcation to a breathing spot observed in the experi- 
ment. 

The oscillations can be regarded as periodic transitions 
between left and right propagating fronts represented by 
the upper and lower branches of the c — d relation in 
Fig. 4(b) [|Io|. The dynamics actually do not follow these 
branches because we have neglected the explicit time de- 
pendence of the u and v fields in the moving frame. Near 
the bifurcation, Vf becomes an active degree of freedom 
responsible for transitions between the fronts JTT[] . The 
present analysis, however, accurately predicts (within 3% 
for ai = 5) the onset of breathing motion and describes 
the dynamics far from the front bifurcation (see Fig. 4(a) 
and 4(c)). 

Beyond the front bifurcation, ao > ao&(e), the upper 
branch in the c — d relation extends to infinity, i.e., ex- 
ists for all distances d. As Fig. 4(c) demonstrates, this 
type of c — d relation results in the rebound of a front 
approaching the boundary along the lower branch (black- 
white front) and the escape to infinity along the upper 
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FIG. 4. The dynamics of fronts near a boundary. The thick 
lines represent the front distance from the boundary, d, and 
front speed, c, computed from the numerical solution of equa- 
tions (|l|) . The thin lines are the solutions to equations (^) and 
(^). (a) Far into the single front regime, a front approaching 
from large x values stops at a fixed distance from the bound- 
ary; ao = —0.2, e = 0.025. (b) Near the front bifurcation 
the c — d relation is multivalued and the front begins oscil- 
lating; ao = —0.1, e = 0.0025. (c) After crossing the front 
bifurcation the upper branch of the c — d relation exists for all 
d and the approaching front rebounds and propagates away 
from the boundary to infinity; ao = —0.01, e = 0.0025. In all 
three cases, a\ — 5 and 5 = 2.0. 

branch (white-black front). This behavior together with 
the circular geometry of the experimental apparatus ex- 
plain the rebound and collapse of black spots observed at 
high flow rates. Possible instabilities to transverse per- 
turbations |?||^] are damped in the present experiment by 
high curvature (small spots) and by the interaction with 
the circular reactor boundary (large spots). The observa- 
tion of traveling rings in the same parameter range as the 
collapsing spots (Fig. 3) provides further evidence for this 
interpretation since crossing the front bifurcation is as- 
sociated with the appearance of traveling waves Q . The 
periodic production of waves near the boundary might be 
due to heterogeneities, which become significant near the 
oscillatory regime at high flow rates, or global coupling 
effects. The latter possibility is not highly probable as 



discussed below. 

We have been able to explain breathing spots in terms 
of dynamic front transitions near a parity breaking front 
bifurcation. The generic nature of this bifurcation sug- 
gests that similar behavior can be expected in other 
systems with fronts. Breathing spots can also arise 
from global coupling, as found by Middya and Luss in 
model equations ||. Indeed, we find that decreasing the 
strength of a global coupling term added to (la) leads 
to a scenario similar to that observed: stationary, oscil- 
lating, and collapsing spots. In our experiment global 
coupling can arise from diffusion of chemicals from the 
gel back into the stirred reservoir. To test whether such 
coupling is significant, we monitor the pH in the reser- 
voir. We find that, despite the oscillations in the gel, the 
pH in the reservoir is time-independent. This suggests 
that the primary mechanism leading to the oscillating 
spots is not global coupling but the one presented in this 
Letter: interactions of the front with the boundary near 
a parity breaking bifurcation. 
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